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ABSTRACT 

Radio pulsar subpulse drifting has been interpreted as rotation of sub-beams (sparks) of pair plasma 
produced by intermittent breakdowns of an inner vacuum gap above the pulsar polar cap. This model 
also predicts strong thermal X-ray emission from the polar cap caused by inflowing particles created 
in spark discharges. We have observed the best-studied drifting pulsar B0943-I-10 with XMM-Newton 
and detected a point source coincident with the radio pulsar position. Its spectrum could be fitted 
with a thermal blackbody model, although a power-law model is also acceptable. The thermal fit gives 
a bolometric luminosity Lboi ~ 5 x 10^®ergs“^ and a surface area A w I0^(T/3 MK)”"^ m^, much 
smaller than the conventional polar cap area, 6 x 10^ m^. Such thermal radiation can be interpreted as 
emitted from footprints of sparks drifting in an inner gap of a height h ^ 0.1-0.2rpc, where rpc is the 
polar cap radius. However, the original vacuum gap model by Ruderman and Sutherland requires some 
modification to reconcile the X-ray and radio data. 

Subject headings: pulsars: individual (PSR B0943-I-I0) - stars: neutron - X-rays: stars 


1. INTRODUCTION 

The mechanism of radio pulsar emission remains a mys¬ 
tery after decades of study (Melrose 2004). It is generally 
agreed that there exists an inner magnetospheric charge 
acceleration region (or inner gap) near the pulsar polar 
cap. Ruderman & Sutherland (1975, hereafter RS75) pro¬ 
posed that the gap is nearly vacuum, based on the as¬ 
sumption that ions cannot be stripped off the surface by 
the electric field component parallel to the magnetic field 
lines. Such an assumption was disputed by later workers 
(e.g., Jones 1986; Kossl et al. 1988), who argued that the 
surface charges can be freely supplied into the magneto¬ 
sphere. The inner accelerator flow is therefore expected 
to be space charge limited (Arons & Scharlemann 1979; 
Harding & Muslimov 1998,2001,2002; Harding et al. 2002). 
Nevertheless, the vacuum gap model is widely used to in¬ 
terpret radio data. Suggestions to amend the binding en¬ 
ergy problem have been made (e.g., Usov & Melrose 1996; 
Xu et al. 1999; Gil & Melikidze 2002), and some radio 
emission models heavily rely upon the existence of such 
an inner vacuum gap (e.g., Qiao & Lin 1998; Melikidze et 
al. 2000). To resolve this controversy, one should look for 
independent observational signatures of the vacuum gaps. 

A distinguishing property of the vacuum gap model is 
its high polar cap heating rate. In this model, the vacuum 
gap breaks down intermittently due to pair production dis¬ 
charges (RS75), and the number of the outflowing parti¬ 
cles is comparable to that of the inflowing particles, which 
results in a high luminosity of polar caps heated up to X- 
ray temperatures (e.g., Zhang et al. 2000). On the other 
hand, the space-charge-limited flow (SCLF) model pre¬ 
dicts a much lower polar cap heating rate because only a 
small fraction of the positrons (for an electron accelerator) 
turn around and bombard the surface (Arons & Scharle¬ 
mann 1979; Zhang & Harding 2000; Harding & Muslimov 
2001,2002). Therefore, measuring the thermal X-ray lumi¬ 


nosity from heated polar caps can provide a clue to verify 
the existence of the vacuum gaps. 

The strongest support for the vacuum gap model comes 
from its ability to explain the regular sub-pulse drifting 
observed in some long-period, old pulsars with the so- 
called “conal” emission beam (Rankin 1986, and references 
therein). The phenomenon is naturally interpreted as the 
E X B drift of polar cap “sparks” circulating around the 
magnetic pole (RS75). One of the best-studied drifting 
pulsars is PSR B0943-kl0 (P = 1.10 s, P = 1.0 x lO^^ 
erg s“^, r = 5.0 Myr, Bp = 4.0 x 10^^ G). Deshpande & 
Rankin (1999, 2001) monitored the drifting pattern and 
revealed the “polar cap map” of this pulsar. They identi¬ 
fied 20 sparks rotating counterclockwise with a period of 
37 rotation periods (i.e., P3 37P, where P is the spin 
period, and P3 is the minimum period for the drifting pat¬ 
tern to repeat itself), a conclusion confirmed by an anal¬ 
ysis of low-frequency radio data (Asgekar & Deshpande 
2001). The result was claimed to be consistent with the 
RS75 model (e.g., Deshpande & Rankin 1999, 2001; Gil & 
Sendyk 2003), although some modifications to the model 
may be needed (Gil et al. 2002; Asseo & Khechinashvili 
2002; Gil et al. 2003). 

According to the Galactic electron density model 
(Gordes & Lazio 2002), the distance to PSR B0943-I-10 
determined from its dispersion measure (DM= 
15.35 pccm“^) is 0.63 ± 0.10 kpc. The expected high 
X-ray luminosity and the small distance make it the best 
target to test the predictions of the vacuum gap polar cap 
heating model. Here we report on observations of PSR 
B0943-I-10 with XMM-Newton, aimed at detecting the 
thermal X-ray emission from the polar cap region. 

2. OBSERVATIONS AND DATA ANALYSIS 

PSR B0943-I-10 was observed twice with the European 
Photon Imaging Gamera (EPIG) aboard XMM-Newton. 
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Fig. 1.— Confidence contours (68%, 90%, and 99%) for the blackbody (left) and power-law (right) model fits to the EPIC-pn spectrum of 
PSR B0943+10. The BB normalization (vertical axis) is the projected emitting area in units of m^, for d = 630 pc. The PL normalization 
is in units of 10“® photons cm“^ s“^ keV“^ at 1 keV. The lines of constant bolometric luminosity (left panel; for a polar cap close to the 
center of stellar disk, in units of 10^® erg s“^) and constant flux (right panel; in units of 10“^® erg cm“^ s“^) are plotted as dashed lines, for 
fixed uh = 4.3 X 10^® cm“^. 


The EPIC-pn (hereafter PN) and two EPIC-MOS (here¬ 
after MOS) detectors were operated in extended full frame 
and full frame modes, respectively. Thin filters were used 
for all the three EPIC detectors. The initial observation 
of 2003 May 7 was interrupted by a strong solar flare. The 
exposure time for this observation, compromised by a high 
background, was 15 ks. Second observation was taken on 
2003 December 2-3, for 20 ks. We processed the data 
using XMM-Newton Science Analysis Software (SAS), v. 
6.0. We removed the intervals of high background and re¬ 
tained only patterns 0-4 and 0-12 for the PN and MOS 
detectors, respectively. We also applied energy filters to 
keep the events only in the energy range 0.5-8.0 keV. The 
resulting good exposure times are 24.3 ks for PN, 32.8 ks 
for MOSl, and 32.9 ks for MOS2. 

We created a composite image using all the data 
(PN+MOS) in the 0.5-8 keV band. Using a modihed 
version of the CIAO wavdetect tooft, we found 22 point 
sources in the 20' x 20' field. We searched for counterparts 
of the field sources in the 2MASS and USNO-Bl.O cata¬ 
logs and found three matches in 2MASS and one match 
in USNO-Bl.O, within 2" of the X-ray source positions. 
(There are no other matches closer than 6".) Assuming 
these are the counterparts of the X-ray sources, we found 
average offsets of the X-ray with respect to optical/NIR 
positions of -l-O'.'l in RA and —0'.'9 in Dec, with an rms er¬ 
ror of about 0'.'5 in each coordinate (which does not include 
the uncertainty of the X-ray source position determina¬ 
tion) . Visual inspection of the image shows a faint source 
at the (boresight-corrected) position a = 09''46™07?4, 

5 = 09°51'54" (J2000.0), with an uncertainty of about 
5" in each coordinate. The main source of the uncertainty 
is a large error in centroiding the very faint X-ray source. 
Since the separation of 2'.'7 between the X-ray source po¬ 
sition and the radio pulsar position (a = 09^46™07?652, 

6 = 09°51'55'.'52) is smaller than the source position un¬ 
certainty, we conclude that the observed source is most 
likely the X-ray counterpart of the radio pulsar. 

^http://cxc.harvard.edu/ciao/threads/wavdetect 


We measured background-subtracted source count rates 
of 2.52 ±0.48, 0.95 ±0.19, and 1.03 ±0.21 counts ks"! for 
the PN, MOSl, and MOS2 detectors, respectively. These 
count rates were measured in the 0.5-8 keV band within 
a 16" radius aperture, which contains « 67% and 71% of 
total point source counts for the PN and MOS detectors, 
respectively. 

We were unable to find a satisfactory way to combine 
the PN and MOS data for spectral analysis of this faint 
source as the responses of these detectors are very differ¬ 
ent. Since the number of source counts in the MOS data 
is very small, we used only the data from the most sen¬ 
sitive PN detector. We extracted the source±background 
spectrum (102 counts) for the combined PN data from a 
16" radius circle and the background from a nearby 32" 
radius source-free region. We binned the spectrum to have 
a minimum of 10 counts per energy bin. Spectral fitting 
was done using the XSPEC package (v. 11.3.0). 

Because of the scarce count statistics, we fixed the ab¬ 
sorbing column to the value corresponding to the pul¬ 
sar’s DM (assuming 10% ionized ISM): nn = 4.3 x 10^° 
cm“^. The spectrum fits equally well with an absorbed 
blackbody (BB) or an absorbed power-law (PL) model 
ixi = 0-9 and 0.7, respectively, for 8 dof). The BB 
model gives a best-fit temperature T = 3.1 MK and a 
projected emitting area of A ~ lO^dlgg m^, correspond¬ 
ing to an effective radius of about ISdeao m of an equiva¬ 
lent emitting disk, where deao = d/(630pc). This ther¬ 
mal radiation could be interpreted as emitted from a 
hot polar cap; in this case the bolometric luminosity is 
Lboi = AaT"^ = x lO^®(i|3Q(cos0)“^ erg s“^, where 

(cos 9) is a time-averaged cosine of the angle between the 
magnetic axis and the line of sight ((cos 6*) = 0.97 for the 
axis orientations inferred by Deshpande & Rankin 2001). 
The errors here and below are determined from the 68% 
confidence contours for joint variation of two fitting pa¬ 
rameters. Although the strongly correlated temperature 
and radius vary in broad ranges, the luminosity is much 
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more restricted, as demonstrated in the left panel of Fig¬ 
ure 1. 

The PL model gives a photon index F = 2.61 q 5 
a flux Fb.5-8keV = d-dtij X 10“^^ erg cm“^ s“^. This 
flux corresponds to an isotropic luminosity Los-gkeV = 
47rd2F^oT-1keV = 2.41^7 x erg s'l. 

The inferred parameters of the X-ray source do not con¬ 
tradict to the hypothesis that it is the counterpart of the 
pulsar (see §3). The most direct way to prove the identi¬ 
fication would be detection of X-ray pulsations. Our tim¬ 
ing analysis of the PN data (time resolution 0.2 s) did 
not reveal any significant signal at the pulsar’s period of 
1.1 s. This result, however, is not restrictive: the radia¬ 
tion should be about 100% pulsed to find pulsations with 
only 102 counts detected in two observations separated by 
about 7 month (of which only about 61 counts are expected 
to come from the source). Thus, given the uncertainty of 
the spectral results and a lack of period detection, we can¬ 
not completely rule out the possibility that we detected 
a field X-ray source very close to the pulsar position. In 
this case, the above flux/luminosity estimates should be 
considered as upper limits. 

3. PHYSICAL IMPLICATIONS 

As described in §2, the quality of the data is too poor 
to discriminate between the thermal (BB from a po¬ 
lar cap) and nonthermal (PL from the magnetosphere) 
spectral fits of the putative X-ray counterpart of PSR 
B0943-I-10. The photon index for the PL fit is somewhat 
larger than, but not inconsistent with, those observed in 
other pulsars (F = 1-2), given its large uncertainty (in 
particular, F decreases if a smaller nn value is assumed). 
The luminosity found from the PL fit allows one to esti¬ 
mate the “X-ray efficiency” of the magnetospheric emis¬ 
sion, rj = Lx/E. It follows from the above results that 
T] = 2.4(-1.2,-Ll.8) X 10-3 and 0.8(-0.6,-f 1.5) x IQ-^ 
in the 0.5-8 and 2-10 keV bands, respectively, where the 
distance uncertainty is included (the latter band is often 
used for characterization of the nonthermal pulsar radia¬ 
tion; e.g., Possenti et al. 2002). Although the inferred ef¬ 
ficiency is rather high (close to the upper limit suggested 
by Possenti et al. 2002), it is within the range of efficien¬ 
cies found recently for other old pulsars (Zavlin & Pavlov 
2004). 

Within the thermal interpretation of the observed ra¬ 
diation, particularly interesting implications follow from 
the comparison with the predictions of polar cap heating 
models, especially the vacuum gap model (RS75). In this 
model, a vacuum gap of height h breaks down due to run¬ 
away electron-positron pair production, forming “sparks” 
along the magnetic field lines. A half of the relativistic par¬ 
ticles produced in the spark discharge stream back toward 
the neutron star and heat the “spark footprints” on the po¬ 
lar cap surface. Since the strong magnetic field suppresses 
the heat conduction across the field lines, and the cooling 
timescales are very short (about a few /rs for B = a, few 
xlO^^G; Gil et al. 2003), only the spark footprints emit 
thermal X-rays while the rest of the polar cap remains cold. 
The BB fit of the observed spectrum allows one to estimate 
the bolometric luminosity, Lboi ~ 5x 10^® ergs s”^, and the 
total X-ray emitting area that is strongly correlated with 
the temperatnre, A « 103(r/3 MK)”^ m^ {A « 300-5000 


m^, T w 2.0-4.2 MK, at 68% confidence). Even the largest 
area allowed by the fit is much smaller than the conven¬ 
tional polar cap area for a dipole magnetic field, ApJP = 
7r(rpJP)^ — 6 X 10^ m^, where TpJP = {2'kR ^/~ 138 
m, for a neutron star radius i? = 10 km. The spark foot¬ 
prints fill a fraction of the polar cap, with a filling factor 
/ = A/Ape ~ O.OI 75 A 3 , where A 3 = A/(103m^), and the 
factor b = Bg/Bp ~ Ap)P/Apc > 1 takes into account that 
the magnetic field is non-dipolar (i.e., the field in the gap, 
Bg, may exceed the dipolar value Bp, and the polar cap 
area may be smaller than Ap)?; Gil et al. 2002, 2003). The 
radio sub-pulse drifting was interpreted as caused by 20 
sparks drifting along the circular outer boundary of the 
polar cap (Deshpande & Rankin 1999,2001). Since some 
sparks may not be seen in the radio (but should be seen 
in more isotropic thermal X-rays), the total number of 
sparks N can exceed 20 (e.g., Gil & Sendyk 2000 suggest 
N ~ 42). A characteristic spark radius, Vgp, can be esti¬ 
mated as (assuming all sparks have the same size) 

rsp - {A/nNy/^ = 4:Al^‘^{N/20)-^/^m. (1) 

Another constraint on the gap and spark parameters 
comes from the comparison of the bolometric and spin- 
down luminosities. The bolometric polar cap luminos¬ 
ity can be estimated as Lboi ~ (l/2)$cp(,jA, where $ ~ 
2ttB gK^/cP is the potential drop across the gap (along the 
magnetic field lines), and ~ Bg/cP is the Goldreich- 
Julian charge density. On the other hand, the total spin- 
down luminosity can be written as E ~ 2<I)i„axCp(i’PAp)? = 
24>maxCPGj-4pc, where ^max — BpAfP/cP = BgApe/cP 
is the maximum possible potential drop along a magnetic 
field line (RS75). This gives Lboi/E ~ (//4)($/$i„ax) = 
(l/2)A^(rsp/rpc)^(h/rpc)^. Using the observational value 
Lho\/E « 5 X 10“^, we obtain 

(h/rpe)(r,p/rpe) - 7 x W-yN/20)-^/^. (2) 

Adopting the assumption that the size of the spark is com¬ 
parable to the gap height (i.e., rsp ~ h/2; Gil & Sendyk 
2000), we obtain from equation (2) 

V?-pc ~0.12(Af/20)-i/% (3) 

This gives h ~ /2Q)-^/'^ m, and rgp ~ 

8 5-%2(7v/20)-i/4 

m. Using this estimate for rsp and 
equation (1), we obtain A ~ 40006-^(A^/20)^/^ m^, which 
is consistent with the X-ray data at reasonable values of b 
and N (see Fig. 1). 

Further constraints can be obtained from the drift¬ 
ing period. According to RS75, for sparks at the ra¬ 
dial distance r± from the pole, it can be estimated as 
P 3 ~ 2'!Tr±/v ~ 2TTr±Bg/Ej_c, where r = |E x Bs|c/i?g 
is the drift velocity, and E± = \dV{r±)/dr±\ is the elec¬ 
tric field component perpendicular to the magnetic field 
lines at the top of the gap. In the case h ^ rpc, U(r_L) 
is approximately uniform throughout the polar cap ex¬ 
cept for near the cap edge, rpc — r± ^ h (Appendix II 
of RS75). This means that the “inner sparks” drift very 
slowly, while Ej_ ~ 4>//i, v ~ 2Trh/P, and P3/P ~ rpc/h 
for the sparks near the cap edge. The measured P3/P ~ 37 
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thus gives hjr^c ~ 0.027, or ft, ~ 4 6“^/^ m and rgp = 
36 ft“^/^(A^/20)~^/^ m (from eq. [2]). Such a spark radius 
substantially exceeds the gap height, and it corresponds 
to the emitting area, A ~ 8 x 10^ft~^ m^, well above the 
range allowed by the X-ray lit. Therefore, the original vac¬ 
uum gap model of RS75 can hardly be reconciled with the 
X-ray data. However, we cannot rule out that a modifi¬ 
cation of the RS75 model could explain both the X-ray 
and radio data in a consistent way. For instance, Gil et al. 
(2003) suggested a modification assuming thermal emis¬ 
sion of charged particles from the polar cap surface which 
can partially screen the gap potential: $ ^ C<I), where 
C < 1 is the screening parameter. In such a modified RS75 
model, Ps/P ~ r-pc/Cft. For ( = 0.17, inferred by Gil et al. 
for PSR B0943-I-10, we obtain ft/rpc ~ 0.16, very close the 
value inferred from the X-ray data (eq. [3]). However, a 
more detailed study of the dynamics of a partially screened 
unsteady polar gap is desirable for an accurate treatment 
of the problem. We should also caution that the BB model 
used for fitting the spectrum gives only crude estimates for 
the area and temperature in the case of strong magnetic 
fields (Zavlin & Pavlov 2004). To conclude, the thermal 
interpretation of the X-ray data is not inconsistent with 
the general idea that the X-ray radiation is emitted from 
spark footprints, but the original vacuum gap model of 
RS75 needs some modification to explain both the X-ray 
and radio data. 

The results of our X-ray observations also constrain 
other models. Motivated by solving the binding energy 
problem, Xu et al. (1999) suggested that PSR 0943-1-10 
and other drifting pulsars might be strange quark stars 
covered by a very thin layer of normal matter. In this case, 
the heat deposited at the polar cap spreads over the en¬ 
tire star surface because of a high thermal conductivity of 
the surface layer (Xu et al. 2001). This results in a nearly 
uniform, low surface temperature (^ 0.1 MK), being in¬ 
consistent with the thermal interpretation of the observed 
X-ray radiation (but of course this hypothesis cannot be 
excluded if the observed radiation is nonthermal). 

Polar cap heating in the SCLF model has been studied 
by Harding & Muslimov (2001,2002). Assuming a dipole 
field configuration near the surface, the timing parameters 
of PSR B0943-I-10 indicate that the gap is controlled by 
inverse Compton scattering (ICS) rather than curvature 
radiation (CR), and that the gap is in the “saturated” 
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regime. Using eqs. (64) and (66) of Harding & Muslimov 
(2002), the predicted polar cap luminosity for this pul¬ 
sar is 8.8(1.7) X 10^^ergs“^ for a resonant (non-resonant) 
ICS gap. Since these luminosities are much lower than 
observed, the SCLF model is inconsistent with the ther¬ 
mal interpretation of the X-ray data, but is consistent with 
the nonthermal interpretation. If future observations show 
that the X-ray radiation is, in fact, predominantly non¬ 
thermal (Lboi ^ 5 X 10^® ergs s“^), then the regular drift¬ 
ing behavior observed in PSR B0943-I-10 and other drifting 
pulsars will have to be interpreted by models other than 
bunched plasma sparks in the inner gap (e.g., Kazbegi et 
al. 1996; Wright 2003). 

4. SUMMARY 

To test the existence of an inner vacuum gap in PSR 
B0943-I-10, we observed the pulsar with XMM-Newton 
and detected a point source at the radio pulsar posi¬ 
tion. The best-fit 0.5-8 keV band isotropic luminosity is 
~ 2.4 X 10^® erg s“^ for the PL fit, while the bolometric 
BB luminosity is about 5 x 10^® erg s“^. Within the ther¬ 
mal interpretation, the X-ray radiation is emitted from a 
heated area much smaller than the conventional polar cap 
area. In the framework of the sparking pulsar model, this 
implies that only spark footprints are heated up to X-ray 
temperatures while the rest of the polar cap remains rel¬ 
atively cold. The thermal model, supplemented with the 
results of observations of subpulse drifting, is consistent 
with the presence of an inner gap of a height ft ~ 0.1- 
0.2 Tpc, and fV > 20 rotating sparks whose footprints fill 
~ 0.04-0.07 of the polar cap area. The thermal model 
disfavors the strange quark star model of PSR 0943-1-10, 
and it is also inconsistent with the SCLF polar cap heating 
model. Both models are, however, allowed if the spectrum 
is non-thermal. 

The nonthermal (magnetospheric) interpretation of the 
observed radiation cannot be ruled out. In this interpreta¬ 
tion, the X-ray efficiency, rj ~ 10“® in the 2-10 keV band, 
is comparable to that of other old pulsars. Much deeper 
observations are required to firmly establish the nature of 
the X-ray radiation from PSR B0943-I-10. 
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